A STUDY OF THE NAVIER-STOKES EQUATIONS WITH THE 
KINEMATIC AND NAVIER BOUNDARY CONDITIONS 



GUI-QIANG CHEN AND ZHONGMIN QIAN 

Abstract. We study the initial-boundary value problem of the Navier-Stokes equations 
for incompressible fluids in a domain in with compact and smooth boundary, sub- 
ject to the kinematic and Navier boundary conditions. We first reformulate the Navier 
boundary condition in terms of the vorticity, which is motivated by the Ifodge theory on 
manifolds with boundary from the viewpoint of differential geometry, and establish basic 
elliptic estimates for vector fields subject to the kinematic and Navier boundary condi- 
tions. Then we develop a spectral theory of the Stokes operator acting on divergence-free 
vector fields on a domain with the kinematic and Navier boundary conditions. Finally, 
we employ the spectral theory and the necessary estimates to construct the Galerkin 
approximate solutions and establish their convergence to global weak solutions, as well 
as local strong solutions, of the initial-boundary problem. Furthermore, we show as a 
corollary that, when the slip length tends to zero, the weak solutions constructed con- 
verge to a solution to the incompressible Navier-Stokes equations subject to the no-slip 
boundary condition for almost all time. The inviscid limit of the strong solutions to the 
unique solutions of the initial-boundary value problem with the slip boundary condition 
for the Euler equations is also established. 



1. Introduction 

We are concerned with solutions of the initial-boundary value problem of the Navier- 
Stokes equations for incompressible fluids in a general domain in with compact and 
smooth boundary, subject to the kinematic boundary condition (i.e. the slip condition) 
and the Navier boundary condition. The incompressible fluid flows are governed by the 
Navier-Stokes equations: 

dtu + u-Vu = fj,Au-Vp, V-u = 0, a; G C M^, (1.1) 

where u represents the Eulerian velocity vector field of the fluid flow, p is a scalar pressure 
function (up to a function of time t) which maintains the incompressibility of the fluid, n 
is the kinematic viscosity, and the density has been renormalized as one in this setting. As 
a system of partial differential equations, u and p are unknown functions, and a solution u 
of (jl.ip determines p uniquely up to a function depending only on t. The initial condition 
for the fluid flow is 

u\t=o = uo{x). (1.2) 
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If C M'^ has a non-empty boundary F = 50, then system (jl.ip must be supplemented 
with boundary conditions on u in order to be well-posed. In fluid dynamics, if the rigid 
surface V is at rest, the kinematic and no-slip conditions are often imposed. The kinematic 
condition means that the normal component of the velocity vanishes, that is, the velocity 
u is tangent to the boundary T: 



while the no-slip condition demands for the coincidence of the tangent component of 
the fluid velocity with that of the boundary T. These two boundary conditions lead to 
the Dirichlet boundary problem associated with the Navier-Stokes equations. There has 
been a large literature for the Navier-Stokes equations subject to the Dirichlet boundary 
condition; see [11 [11 EqI |24l ESI [Ml HQI ISl |44] and the references cited therein. The 
fundamental problem of the global (in time) existence and uniqueness of a strong solution 
remains open; however, the Dirichlet boundary problem of the Navier-Stokes equations is 
well-posed at least for a small time, or for small data globally in time. 

However, the usual no-slip assumption does not always match with the experimental 
results. Navier [34] first proposed the slip-with-friction boundary condition, that is, the 
Navier boundary condition: The tangent part of the velocity u is proportional to that 
of the normal vector field of the stress tensor with proportional constant C > Oj which is 
called the slip length (see (j2.5p below). In the recent years, the Navier boundary condition 
has been received much attention, especially when fluids with larger Reynolds number or 
fluids past a rigid surface with considerable speeds for which the curvature effect becomes 
apparent (cf. [HI [151 113 [13 [231 [Hj and the references cited therein). Such boundary 
conditions can be induced by effects of free capillary boundaries, a perforated boundary, 
or an exterior electric field (cf. [B [H [3 [H [371 [39]). In particular, this friction boundary 
condition was rigorously justified as the effective boundary condition for flows over rough 
boundary; see y^itlTj. Thus, it becomes important to analyze solutions to the equations 
for such fluids subject to the Navier boundary condition. 

The rigorous mathematical analysis of the Navier-Stokes equations with the kinematic 
and Navier boundary conditions may date back the work by Solonnikov-Scadilov [12] for 
the stationary linearized Navier-Stokes system under the boundary condition that the 
tangent part of the normal vector field of the stress tensor is zero. The existence of weak 
solutions and regularity for the stationary Navier-Stokes equations with the kinematic 
and Navier boundary conditions was only recently obtained by Beirao da Veigt [6j for the 
half-space. In a two-dimensional, simply connected, bounded domain, the well-posedness 
problem has been rigorously established by Yodovich [l6]. See also Clopeau, Mikelic, and 
Robert [9] and Lopes Filho, Nussenzveig Lopes and Planas [28] for the vanishing viscosity 
limit, and Mucha [33j under some geometrical constraints on the shape of the domain. 
These two-dimensional results are based on the fact that the vorticity is scalar and satisfies 
the maximum principle. However, in the three-dimensional case, the standard maximum 
principle for the vorticity fails, so that the techniques employed in the two-dimensional case 
can not be directly extended to this case. Furthermore, the Navier boundary condition 
causes additional difficulties in developing apriori estimates which require to be compatible 
with the nonlinear convection term. The main purpose of this paper is to develop a general 
approach to establish the well-posedness, the no-slip limit, as well as the inviscid limit for 
the initial-boundary value problem for the Navier-Stokes equations in a general domain 
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C M'^, subject to the Navier boundary condition, together with the kinematic condition 
(jl.3p and the initial condition (jl.2p . 

By careful local computations in Section 2, we first reformulate the Navier boundary 
condition in terms of the vorticity u; = V x u (see Proposition 12.11 below). This is moti- 
vated by the Hodge theory on manifolds with boundary, since the kinematic and vorticity 
conditions are the natural boundary conditions for the Hodge theory from the viewpoint 
of differential geometry. Indeed, the Navier boundary condition under the kinematic con- 
dition (jl.3p is equivalent to the condition that the tangent portion of the vorticity to is of 
the following form: 



where vr = (vTjj) is the curvature of the boundary F, that is, the second fundamental form. 
The form vr is identified with the self-adjoint operator which sends a tangent vector 
on the boundary surface F to the tangent vector (^^- T^iji,-' ,Ylj '^2j(,-')j the operator * is the 
Hodge star operator which rotates (towards the interior of the domain Q) a tangent vector 
(C^) by 90*^ degree. The Navier boundary condition in form (II. 4|) has appealing physical 
interpretation: The vorticity on the boundary F is created mainly from the slip of the fluid 
and the curvature of the boundary F (see [23] for more information and further references 
on the slip length for different fluid media). In particular, the effect of the curvature 
becomes significant when the curvature of the boundary F becomes large (comparable to 
the reciprocal of the slip length). 

In Section 3, we establish some basic elliptic estimates for vector fields subject to the 
kinematic and Navier boundary conditions. In particular, we establish some L^-estimates 
which are uniform in the slip length by a careful analysis of several boundary integrals. 

In Section 4, we develop a spectral theory of the Stokes operator acting on divergence- 
free vector fields in a general domain $7 subject to the kinematic and Navier boundary 
conditions. We establish several fundamental estimates for the symmetric form defined 
by the Stokes operator. Besides the difficulties caused by the divergence-free condition on 
the vector fields, the Navier boundary condition causes additional difficulties in developing 
apriori estimates for the Galerkin approximations to solutions of the Navier-Stokes equa- 
tions. To overcome these difficulties, we establish an estimate for the third derivatives of 
the vector fields satisfying the Navier boundary condition (Theorem 14.41 and Corollarv l4.2p 
and a uniform gradient estimate for the Galerkin approximations (Theorem 14. 5p . Then, 
in Section 5, we employ the spectral theory and all the estimates established in Sections 
3-4 to construct the Galerkin approximate solutions and establish the global existence of 
weak solutions and the local existence of strong solutions for the initial-boundary problem 
(jl.ip ~ (jl.4p . Furthermore, we show as a corollary that, for any weak solution M^(t, x) cor- 
responding to the slip length ( to problem (jl.ip - (jl.4p constructed in Theorem 15.11 when 

— > 0, there exists a subsequence (still denoted) n^(t, x) converging to u{t,x) such that 
u{t,x) is a solution to (jl.ip subject to the no-slip condition for almost all time t. On the 
other hand, when ( — > oo, there also exists a subsequence (still denoted) ti^(t, x) converg- 
ing to u{t, x) such that u{t, x) is a solution to (jl.lj) subject to the complete slip boundary 
condition: 




(1.4) 
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in the weak sense. Such a nonhomogeneous vorticity boundary problem has been carefully 
investigated in [7j. 

Finally, in Section 6, we study the inviscid limit and establish the L^-convergence of 
the strong solutions of problem (jl.ip - ()1.4p to the unique smooth solution of the initial- 
boundary value problem with the slip boundary condition for the Euler equations for 
incompressible fluid flows. 

2. The Navier boundary condition 

In this section we reformulate the Navier boundary condition in terms of the vorticity 
for every slip length C > introduce several notions and notations which are used 
throughout the paper. 

For simplicity, we use the conventional notation that the repeated indices in a formula 
are understood to be summed up from 1 to 3 unless confusion may occur. Furthermore, 
we use a universal constant C > that is independent of the slip length ^ > 0, and a 
universal constant M > that may depend on ^ among others, which may be different at 
each occurrence. 

We use the same notation for both scalar functions and vector fields in the L^-space 
and Sobolev spaces W^'P{^) {H^{^) ii p = 2). Denote ||T||p as the L^-norm of the length 
|T| of T on O with respect to the Lebesgue measure, and ||T||^p(p) as the L^-norm of the 
vector field T on the boundary F with respect to the induced surface-area measure on F. 
That is, 

IITIIp = ( 1^ \T{xW dxf, ||T|Up(r) = ( \T{xW dn\x)f, 

where dx is the usual Lebesgue measure on and is the two-dimensional Hausdorff 
measure (i.e. the surface area measure) on F. From now on, dx and dTi.'^{x) in the integrals 
will be suppressed, unless confusion may arise. For a vector field T on 

||T||h.... = (^ / N'T\^)K 

where V^T is the j-th derivative of T. See [21 [13] for the details. 

For the bounded domain $7 C with a boundary F = dQ that is smooth, compact, 
and oriented, unless otherwise specified, we carry out local computations on the boundary 
in a moving frame compatible to F. More precisely, if u is the unit normal to F pointing 
outwards with respect to $7, by a moving frame we mean any local orthonormal basis 
(ei, 62, 63) of the tangent space TQ such that 63 = u when restricted to F. If n = ^^=1 u^ej 
is a vector field on Q, then, restricted to the boundary surface F, u" = Y^j=i 2 ^"'^j (resp. 
u-L = u^u) denotes its tangent part (resp. normal part). The Christoffel symbols F-^ 
are determined by the directional derivatives VjCj = rfj-e/cj where V, is the directional 
derivative in the direction Cj. 

The tensor {T^ij)i<i.j<2, where iTij = — Ff^ for i,j = 1,2, is a symmetric tensor on F, 
which is the second fundamental form, denoted by vr. That is, 

7r(u",f") := TTijU^v^ for any u",?;" G TV. 
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We say that vr is bounded above (resp. below) by a constant A if two eigenvalues (which 
are functions on T) are bounded above (resp. below) by A. We will also identify tt with 
the linear transformation (hij): If = 2 ^"'^j is tangent to T, then 

7r(ull) := 7r(u")-'ej = TTijU^ej, 
i=i,2 j=l,2 

with (7r(Mll), w") = 7r(iill , wll), and H = X]j=i2^ii mean curvature. We refer to 

[35\ [38] for further facts and notations in differential geometry used in this paper. 

The boundary surface F C has a natural induced metric and hence a natural notion of 
directional derivatives, the Levi-Civita connection, denoted by V'". The following formulas 
will be useful in treating with integrals on the boundary T. Let u £ H^{i^) be a vector 
field on fi. Then, on T, 

{u-Vu,v) = —it{v}\v}^) — H\u^\'^ + {u,v)V ■ u 

+2(nll,V^(n,zy)) - • {{u,u)J), (2.1) 

and 

-du{\u\^) = (n X (V X n), I/) + (n • Vn, v). (2.2) 

The first formula (j2.ip may be verified by means of the moving frame method (see [7J for 
the details). The second follows from the vector identity: 

lv|np = n X (V X -u) +n • Vii. (2.3) 

Let / be a scalar function on Vt. Then, as a special case of ()2.2p . 

5.(|V/|2) = -27r(V/ll,V/ll)-2F|a,/|Va,/A/ 

+4(v/ii, vr(a,/)) - 2vr . (a./yr/ii). (2.4) 

The connecting condition over an interface F of a fluid is expressed as the Navier bound- 
ary condition; see Einzel-Panzer-Liu [llj for its physical interpretation. This condition 
may be written in a moving frame compatible to F as follows: 

u'' = -C (Vgn'^ + Vku^) on F for A; = 1, 2, (2.5) 



where C is the slip length that is a positive scalar function on F depending only on the 
nature of the fluid and the material of the rigid boundary. In order to write down (12. 5p 
in a global form in terms of the vorticity and the curvature of F, we recall that the Hodge 
operator * sends a vector field {v^,v'^) on the surface F to *{v^,v'^) := {—v'^,v^). The 
effect of the Hodge operator * is to rotate a vector on F by 90" degree with respect to 
the normal vector pointing the interior of ^l. The Hodge operator * is independent of the 
choice of a moving frame on F and may be defined via the identity: 

{w X {*u^,u) = {u\w^^) on F (2.6) 

for any vector fields u and w. 

Proposition 2.1. Let u G C^{Q,). Then u satisfies the Navier boundary condition on F 
if and only if 

(Vxm)II ^ = -i(*iill) -2(* V^(m,i/)) +2(*7r(nll)). (2.7) 
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In terms of the components in a moving frame compatible to the boundary, condition (|2.7p 
takes the following forms: 

(V X u)^ = ^u^ + 2V2(n, I/) - 2 ^ Trj2u\ 

(V X uf = -^u^ -2Vi{u,u) + 2 TTjiuK 

In particular, if u satisfies the kinematic condition (|1.3p . then the Navier condition ()2.5p 
is equivalent to (jl.4p . 

Proof. It suffices to show the results in a moving frame compatible to the boundary surface 
r. Then, for k = 1,2, 

VkU^ = ek{u^) + ^ rljU^ = ek{u, J^) - ^ vr^fcU^. 

i=i,2 i=l,2 

Therefore, along the boundary surface T, 

V'iU^ + VkU^ = e-:ika^°- + 2ek{u, i/) - 2 ^ iijkU^ for any /c = 1, 2, 

i=i,2 

so that 

w^lr = -TEsfcaw'' - 2e3kaek{u, I') + 2ezka ^ T^jkU^ for a = 1, 2, 
^ i=i,2 

where £ijk is the Kronecker symbols. This completes the proof. □ 

Definition 2.1. Let C > be a constant. Then a vector field n on is said to satisfy the 
Navier's (^-condition if (jl.4p holds, which is equivalent to 

{V X uf = ^u'^ -2^ ■Kj2u\ (V X n)^ = -^n^ + 2 ^ TTji-u^' (2.8) 
^ i=i,2 ^ j=l,2 

in a moving frame compatible to the boundary surface T. 

In this paper we study the initial-boundary problem (jl.ip - (jl.4p for the Navier-Stokes 
equations with fixed constants fi > and C > 0. 

3. Elliptic estimates for vector fields 

The fundamental estimates in the standard elliptic theory state that, for any function / 
on Q subject to a certain boundary condition (Dirichlet or Neumann), ||/||f/2 is dominated 
by the L^-norm of its Laplacian together with its L^-norm: 

11/11^. <C7(||A/||2 + ||/||2) 

for some constant C depending only on Q. A correct boundary condition here plays 
an essential role, and the previous estimate can not be true without a proper boundary 
condition. A version of elliptic estimates for vector fields has been established in p] (also 
see |32j ) for vector fields satisfying the Dirichlet or Neumann condition. If n is a vector 
field in H'^{Vl) such that 'U"'"|p^ = and u\y^ = for F = Ti U then 

\\u\\m{n) < C'diV X u\\2 + ||V • u\\2 + \\u\\2) , (3.1) 
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which is a special case of a general result in [^. 

For our problem, we need to develop the L^-estimates for divergence-free vector fields 
that satisfy the kinematic condition (II. 3p and Navier's ^-condition (jl.4p in a general do- 
main Q. To our knowledge, these estimates are not covered in the previous literature, 
although they can be considered as a part of the standard elliptic theory. 

We begin with an elementary lemma which implies (jS.ip and may be verified by means 
of integration by parts. 

Lemma 3.1. Let u G H'^{Q) he a vector field on 17. Then 

{^u,u) = -\\Vu\\l + ]^ j du{\u\^) (3.2) 
and 

|Viip = ||V X u\\l + ||V • u\\l - j (V • u) (n, v) + j {u ■ Vn, i/). (3.3) 
In particular, if u-^\^ = 0, then 

||Vu||2 = ||V X -ull^ + ||V • n||2 - y 7r(n,n). (3.4) 
Lemma 3.2. If g is a smooth function on Q (up to the boundary T), then 

||V25||2 = ||A<7||i-yvr((V<7)",(V5)ll)-yF|a,<7|' + 2y(Vr5,Vr(5,5)). (3.5) 

This elementary fact can be proved by using integration by parts and the Bochner 
identity: 

\V^g\^ = ^A\Vg\^-{VAg,Vg). (3.6) 

Now we are in a position to prove our first main estimate, an elliptic estimate, for vector 
fields satisfying (fOll - ffLl]) . 

Theorem 3.1. There exists a constant C depending only on such that 

+ < C {\\Au\\l + ||n||^,) (3.7) 

for any vector field u satisfying p.3p ~ (!1.4p . 

Proof. Under an orthonormal frame, we have 

3 



k=l-^^ i,j=l,2 
3 

- / FV|a,nY + 2 / (V^n^V^(a,n'=)) 

3 

\\Aug - [ 7r((Vn^')ll,(V^x'=)ll)- / FV(Vn^z.)^ 
+2 J (V^^z^V^(V^x^^y)), 



(3.8) 
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where the second equahty follows from (j3.5p applying to g = u^. The second and third 
boundary integrals can be dominated by J-p |Vup. Thus, we have to handle the last 
boundary integral, where we use the Navier boundary condition (II. 4p . Working in a frame 
compatible to T, since w"'"|p = so that V^u^ = 0, then 

On r, and are the tangent components of u, and = 0, 

d^u'' = e3(n^) = Vsn" - n^T^,, 

and VfcU^ = — X]i=i 2 '^^^ki, k = 1,2. For w = V x u, Vsu'^ — Vjtu^ = s^kj^'' so that, for 
A; = 1,2, 

S^^i'' = £3kjUJ^ + Vfcn^ - u^T^j = eskjUJ^ - ^ u'ttm - n^T^j. (3.9) 

j=l,2 i=l,2 j=l,2 

According to the Navier 's (^-condition ()1.4p (also see the proof of Proposition 12. ip : 

^ b=l,2 

Substitution it into ()3.9p yields 

duu'' = -^eSajeSkjU" + '^£3kj£3aj Y ^baU^ ~ ^ U^ki - ^ U^T^j 

6=1,2 i=l,2 j=l,2 

= -7^' + E ^^'^^^ - E ^'^3, • (3.10) 

i=l,2 j=l,2 



It follows that 



so that 



^ i=l,2 j=l,2 

k=l,2 

= -^iv^^p + E (v^^', v^( E ^^^-') - E E ^^Ts,)) 

^ A:=l,2 i=l,2 k=l,2 j=l,2 

<_i|Vr^x|2 + c(|Vrn|2 + |^x|2). (3.11) 



Therefore, we have 



/<_!/■ |vr^,|2 + C / (|Vn|2 + |^,|2), 
C Jr Jr 

where C is a constant depending only on fi. Combining this inequality with (j3.8p yields 

+ ^ [ |V^n|2 < ||Au||| + C / (|Vnp + 
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Then the conclusion fohows from the Sobolev embedding: 

/ |Vup < ellV^nlli + -||Vn||| (3.12) 

for some constant C = C(r2), independent of e > 0, since the boundary T has bounded 
geometry. □ 

As a consequence, we have the following elliptic estimate. 
Corollary 3.1. There exists a positive constant C depending only on VI such that 

M]i2 + ^||V''n|||2(p) < C||(V X (V X n), V X u, u)\\l (3.13) 

for any vector field u G H'^{^) satisfying (|1.3p - (|1.4p . 

Therefore, for any divergence-free vector field u on 0, satisfying (jl.3p ~ (jl.4p . 

C\\{Au, V xu, u)g < \\u\\jj2 < C~^\\{Au, V xu, u)\\l (3.14) 
for some constant C > depending only on the domain Q, but independent of C > 0. 

4. The Stokes operator with the Navier Boundary Condition 

In this section, we develop a theory of the Stokes operator acting on divergence-free vec- 
tor fields in a general domain $7 subject to the kinematic and Navier boundary conditions 

Note that the divergence operator V- defined for smooth vector fields with compact 
supports in O is closable in L^(Q). The kernel, ker(V-), is a closed subspace of L'^{Q), 
denoted by K2{^). Any vector field u S i^2(^) H H^{^) is divergence-free: V • u = 0, and 
satisfies the kinematic condition (jl.3p . The orthogonal complement of K2{0,) is a closed 
subspace of L^(0), denoted by G2{^), and the decomposition 

is called the Helmholtz decomposition. Any element in G2{^) can be identified with the 
gradient of a scalar function, that is, 

G2m = {vp E L\n) : p e Li^m- 

Let 

Poo : L\n)^K2in) 

be the projection from L^(J7) onto K2{^). The following fact is easy but important. 

Lemma 4.1. Let u G H^{Q), and let u = -Poo(^) + ^Qu be the Helmholtz decomposition 
of u. Then 

V X Poo(m) = V X u, V-Poo(m) = 0, Poo(^^)^|^ = 0. 
Proposition 4.1. Let u £ H^{U). Then 

||VPoo(n)||2 = l|V X ug - j^TT{P^{u),Pooiu)), (4.1) 

and 

||VPoo(n)||2<C||(Vxu, n)||2 (4.2) 
for some constant G > depending only on ^l. 
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Proof. By Lemma |4. 11 we have V x Poo{u) = V x u. Following the elliptic estimate (j3.ip . 

IIVP00WII2 < C\\{V X Poo(n), PooM)||2 < C\\{V X U,U)\\2, 

which gives (j4.2p . Then ()4.ip follows from integration by parts and (j4.2p . □ 

The Stokes operator S can be defined to be the composition S = Poo ° ^ with domain 
We often restrict the Stokes operator on the Hilbert space K2{0,), hence with 
domain H'^{n) n K2{il), but we will use the same notation S if no confusion may arise. 

4.1. The Stokes operator vi^ith the Navier boundary condition ()1.4p . Let C > 

be a constant, and let Dq ,^{S) be the space of all vector fields u G i^2(^) H C°°{Q) (so 
that V-n = and 'u-'"|p = 0) which satisfy the Navier 's (^-condition ()1.4p . Then Dq^(^{S) is 
dense in Er2(i7), and {S, Do,c(5')) is a densely defined linear operator on the Hilbert space 

Lemma 4.2. Let u G Dq^(^(S), and let Au = S{u) + Vp be the Helmholtz decomposition of 
Au. Then p is the unique solution (up to a constant) of the Neumann boundary problem: 

Ap = V-{Au), (9^p|p = iv^-ti-2V^-7r(n). (4.3) 

Proof. Since Au = Su + Vp, then, by taking the divergence and considering the normal 
components, we can easily see that p satisfies the Poisson equation: 

Ap = V-(An), dup\^ = {Au,u). (4.4) 

Since V • n = 0, An = —V x (V x u) so that 

(An, I/) = -(V X (V X u),i^) = -V^ X (V X n)ll, 

and ()4.3p follows from the Navier 's ^-condition (jl.4p . □ 

Remark 4.1. For any p > 2, there exists a constant C{p) > depending only on Q (e.g., 
C(2) = 1) such that 

\\Vp\\p < C\\Au\\p. 

Therefore, 

\\Siu)\\p < C(p)||An||p for any n G Do,^{S). 
Of course, ||5(n)||2 = H-PooA^^lb < ||An||2, if n G H^{n). 
Theorem 4.1. Consider the bilinear form , Dqj^(S)) on K2{i^): 



£{u,w) = — {Su,w) for any u,w E Dq(^{S). (4.5) 
Jn 

Then 

(i) The bilinear form , Dq^^{S)) on the Hilbert space K2{^1) is densely definite, sym- 
metric, and 

£{u,w) = I {Vu,Vw) + - I {u,w) — I ■7r{u,w) for any u,w G Do^(^{S). (4.6) 
Jn C Jr Jr 

(ii) For any e G (0, 1), there exists a constant C(e, $7) such that 

£{u,u)>{l-e)\\Vu\\l-C{e,n)\\u\\l for any u e Dq^^{S). (4.7) 
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(iii) {£, Dq^(^{S)) is closable on K2{0,), its closure is denoted by {£, D(^{£)). Identity 
(|4.6p remains true for any u,w £ Dq{£). 

(iv) If IT < ^, then 

£{u,u)>\\Vu\\l for anyue D(;{£). (4.8) 

(v) Dq{£) = K2{^) n H^{Q.) which is thus independent of C, and hence denoted by 
D{£). 

Proof. Let u,w £ Dq^(^{S) and write Au = S{u) + Vp. Since V ■ u = 0, then 

S{u) = -V X (V X u) - Vp, (4.9) 

where p solves the Neumann problem ()4.3p . Taking inner product on both sides of (14. 9p 
with w and integration by parts on J7 yields 

£{u,w) = {V X {V X u),w) + {Vp,w) 



(V X u,V X w) — 2 I 7r(n, w) + ^ I {u, w), 



{V X u,V X w) + I ((V X n)'l X w, v) 

1 

c 

where we have used (jl.3p - (jl.4p so that 

((V X u)ll X v) = ^(u, w) — 27r(u, w). (4-10) 

Therefore, {u,w) £{u,w) is symmetric and bilinear. Since u-^\^ = w-^\y- = 0, then 

/ {'Vu,Vw) = / {V X u,V X w) - / tt{u,w), 
Jq Jn Jr 

and hence 



£{u,w) = [ {Vu,Vw) + ^ [ {u,w) - [ Tr{u,w). (4.11) 
Jn C Jr Jr 



If vr < ^, then 



£{u,u) > \\Vuf 



12 • 

Let Ai be a upper bound of the second fundamental form vr, i.e., vr < Ai, then 



£{u,u) > \\Vu\\l - Ai / |up > (1 - e) ||Vu||2 ||n||2 

Jr e 

for some C = (7(17) > 0, where we have used the trace imbedding inequality: 

|up < e||Vu||| + ^||n|||. (4.12) 

Next, we show that {£, Z)o,c('S')) is closable on K2{0,). Indeed, if Un G Dq^(^{S) such that 
£{u n — Ujn,Un ~ ^m) — ^ and — timlb — ^ Oj then, since 

we have 

||V(n„ - Ura)\\l 0. 



u 



n~*oo 

SO that 
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That is, {un} is a Cauchy sequence in H^{i^), and hence there exists a unique u £ H^{^}) 
such that 

ll^^n - u\\l + ||V('U„ - u)\\l 0. 

It follows by the Sobolev imbedding that 

lim / {un,Un) = / \uf, lim / 7r(u„,M„)= / 7r(n, 

lim £{un,Un) = [ |Vn|^ + t / - / 7r(u, n), 

and u belongs to the closure of {£, Dq (^{S)). 

Finally, we prove that Dq{£) = K2{^) H H^{fl), which is not surprising, since the 
Navier's ^-condition (II. 4p that has to be satisfied for any u G Dq^(^{S) will be "forgot" 
when passing to the limit in H^(Q) (in which the boundary values of the first derivative 
can not be retained). Therefore, Dq^(^{S) is dense in K2{^) H H^{Q) in the ff^-norm. 

To see this, consider the case that Q = {{xi, X2, x^) : X3 > 0}, and u = (n^,n^,n^) G 
Do^({S) n Cq (M^) such that V ■ u = and u'^{xi,X2, 0) = 0. Then, for every e > 0, choose 

('U^,u^,it^) = (n\n^u^) +X3A'||^3[<g}(/ii,/i2,/i3). 

Then 

{Vxue)^\^ = dx^u^l^ = dx^u'^{xi,X2,0) + h2{xi,X2), 
{Vxuef\^ = - d^^u^l^ = -dx3U^{xi,X2,0) - hi{xi,X2), 
V -Ue = '^{|a;3|<£} (^^sV ■ h + h^) . 

To match the Navier's C-condition, we set 

1 

c'' 

h3ixi,X2) = -X3 (9^1 /ii + 5^3/12) . 
Hence, G D(^{£) and tig ^ n in H^, which concludes the proof. □ 

Corollary 4.1. {£,D{£)) is a densely defined, bounded below, and closed symmetric form 
on the Hilbert space K2{i^). Moreover, 

£{u,u) = ||Vti||2 + ^||^^|Il2(p) — J 7r{u,u) for any u e D(^{£), 

and there exists M{e, Q > such that 

\\Vu\\l < {£ + AI){u,u) < (l + e)||Vu||^ +M(e,C)||n||i for all u £ D^{£). (4.13) 
Proof. Suppose that vr > —Co for some Co > 0. Then 

£{u,u) < \\Vu\\l + + C'o)lhlli2(r)- 

The Sobolev imbedding yields that, for every e G (0, 1), there exists Co > such that 

(i + Co)||u||22(p) < e\\Vu\\l + M(e, C)||n||2 

so that (|iJ3]) follows. □ 



hj{xi,X2) = -U^{xi,X2,()) - d^-iU^{xi,X2,Q), J = 1,2, 
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Definition 4.1. Let ^ > 0. Then the unique self-adjoint operator on K2{^1) associated 
with the closed symmetric form {£, Z)^(<S)) is denoted again by S, with its domain D(^{S), 
called the Stokes operator with Navier's ("-condition, or simply the Stokes operator if no 
confusion may arise. 

According to Definition 14. H (S, D(^{S)) is the unique self-adjoint operator on K2{^) 
such that 



£{u,w) = — {Su,w) for any u G D(^{S), w € D(^{S), 
Jn 

and 

Moreover, if n G D(^{S), then u € H^{il.) with V -u = and n"'-|p = 0. In particular, there 
exists A > such that —S + AI is positive definite (when tt < ^, A = 0). 

To end this section, we establish an L^-estimate for the total derivative of S{u). 

Lemma 4.3. Let Pu be the unique solution (up to a constant) of the Neumann problem: 
Apu = 0, d^pu\^ = ■ u - ■ Tr{u) (4.14) 

for u € H'^{Q) satisfying U"'-|p = 0. Then, for every s > 0, there exists a constant M{e, C) 
depending only on e, Q, and the domain such that 

W^Puh < e||V X V X n||2 + M(e,C)||n||2. (4.15) 

Proof. By integration by parts, one obtains 

IIVpjl = -/p.Ap„WpAp„ 

Jo. Jt 

= ^^p„V^-n-2^p„V^-7r(n) 

< IKIlL2(r)^^(e||V^n||2 + \\u\\h^) 

< \\^Pu\\l^{t) (e||V\||2 + \\u\\Hi) 

for some constant M > which may depend on C, and e. Then (j4.15p follows from (|3.ip . □ 
Proposition 4.2. Let u G Dq (^{S), lo = V x u, and ip = V x uj = —An. Then 

||5(n)||^i <M(||VxV'||2 + ||(V',«)l|2), (4.16) 
where M(e,C) > depend only on e, Q, and il. 
Proof. Since V • S{u) = and S'(n)^ |p = 0, according to p.ip . 

\\S{u)\\l, <C\\{S{u),V X S{u))\\l 

where C > depends only on 0. 

Let S{u) = — — Vp, where p solves (I4.14p . Then V x S{u) = —V x ^ so that 

\\S{u)\\l.=C\\{Vxi^,S{u))\\l 

Together with ()4.15p . ()4.16p follows immediately. □ 
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4.2. Spectral theory of the Stokes operator subject to the kinematic and Navier 
boundary conditions (jl.3|) — (|1.4p . To establish other important properties of the Stokes 
operator (5, we study the boundary problem of the Stokes equation: 

Xu + Au-Vp = f, V-n = (4.17) 

subject to the boundary conditions (fL3]) - (fLl]) . where / G K2i^) n C°°(0) and A e M is 
a constant. 

Taking the divergence to both sides of equation (j4.17|) yields that the scalar function p 
satisfies the Neumann problem of the Laplace equation: 

Ap = 0, dup\-^ = ^V^ -u-lV^ -Triu). (4.18) 

Let A > be the constant such that —S + AI > 0. Then, for A > A, let R\ denote the 
resolvent, i.e., R\ = {XI — S)~^, which is a bounded linear operator on K2{^1). 

Theorem 4.2. For any A > A, Rx is a compact operator on K2{0,). 

Proof. Let / G K2{n) and u = Rxf. Then u G Dt;{S) and 

{XI-S)u = f. 

Suppose in addition that u £ Dq^i^{S), so that V • u = 0, and u satisfies (11.3j) - (ll.4p . Let 
Au = Su + Vp. Then 

Xu- Au + Vp = f, (4.19) 
and p solves the Neumann problem (j4.18p . Hence, for every e > 0, 

\\Vp\\2<e\\V^u\\2+M{eX)\\u\\H^ 
for some constant M(e, ("). It is easy to devise the energy estimate for u. Indeed, since 

X{u,u) - {u,Su) = {u,f), 

then 

XWuWl - [ {u,Su) = [ {u,f) < WuUfh- 
Jn Jn 

Furthermore, since -5 + A > 0, i.e., - /^(u, Su) > \\Vu\\l - A|| u\\l2 , we have 

(A-A)||n||2 + ||Vu||i< ||n||2||/||2. 

Therefore, we have 

^^^IHI^ + l|Vn||i<^||/||i. (4.20) 

Next, we estimate the second-order derivative of u. Since u satisfies (jl.3p ~ (|1.4p . according 
to the elliptic estimate (Theorem 3.1): 

+ ^\N^^^\\h{^) < C {\\Au\\l + \\u\\]j.) , (4.21) 



then 



2, 



< C{\\Vp\\l + \\fg + \\u\\l.) 

< e||V2n||i + M(e,C)(||n||^i + ||/||^) 
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It follows that 

< M(e,C,A) (||n||^, + ||/||2) < M(e, C, A, A)||/||i, 

so that 

<M(C,A,A)||/||i (4.22) 
for some constant M > depending only on A, and A. Therefore, Rx is compact. □ 

Theorem 4.3. The spectrum of the Stokes operator (S, D(^{S)) with Navier's (^-condition 
(jl.4|) is discrete and belongs to (— cxd,A] for some constant A = A(il,(^). The eigenvalues 
Xj < A can be ordered as 

A > Ao > Ai > • • • > A„ > • • • , A„ i -oo. 

Moreover, there are eigenf unctions an-' San = XnO-n, where an G Dq^i^{S) (so that On 
satisfy (jl.3p - ()1.4p ) such that {a„, : n > 0} is a complete orthonormal basis of K2{^). In 
particular, when vr < ^, A = 0, i.e., Xj < for j = 0, 1, 2, • • • . 

Proof. The standard spectral theory for self-adjoint operators yields that the spectrum 
of (S,D(^{S)) belongs to (— oo,A] and is discrete, and there exists an orthonormal basis 
{an '■ n > 0}, where each a„ G and San = ^nO^n for the corresponding eigenvalues 

A„. The standard elliptic theory then implies that o„ G Dq^(^{S). □ 

4.3. Several facts about the Navier's ^-condition (jl.4p . In what follows, we will 
assume that {a„ : n > 0} is the orthonormal basis of X2(0) constructed in Theorem 14. 3[ 
Let is an integer, which may be infinity. Let Xtv be the Hilbert space spanned by 
{ofc : k < N}, and let P/v : L^(r2) — > be the projection. That is, for every u G L^{Q,), 

N 

PNU = \2^k / {ak,u). (4.23) 

Of course, PooU = SfcLo^^ In^'^kju) is the projection from L^(i7) onto K2{^). If < oo, 
Pn{u) G DoxiS) for any u G L'^{n). 

Proposition 4.3. Let u G UN£fqX]\f the vector space spanned by elements in Xjsf where 
N runs over all natural numbers, to = V x u, and = V x u) = —Au. Then 

(V',z^)lr = -Jv^-?x + 2V^-7r(?x), (4.24) 

and 

(V X V)" p = ^( * •5(n)) - 2( * tt{S{u))) . (4.25) 

Proof. Let u G Xjv for some N. By definition, S{u) = Au — Vp so that tp = —S{u) — Vp. 
bmce V • S{u) = and 5(u)-L|p = 0, then p is the solution of the Neumann problem: 

Ap = 0, a^^plp = • n- 2V^ •7r(n), (4.26) 

and {ip, — ~ dup\r, which yields (j4.24p . To see the tangent component of V x "0, we 
note that V x ip = —V x S{u). Since u = Ylk=o^k fni^-kju), then 

N „ iv „ 

S{u) = '^S{ak) / (afc,n) = - ^ AfcCfc / {ak,u) . 
k=o k=o 
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Therefore, 

N 



V X ^ = - V X S{u) = V] Afc (V X ofc) / {ak,u), 



k=0 

and it follows that 

N 



(V X ^ = -^Afc(V X Ofc)" / (afc,u) 
fc=o 

= Afc( -(*afc) - 2( * 7r(aA;)) ) / {ak,u) 

1 ^ r ^ r 



^ k=o A;=0 
= ^{*S{u))-2{*7r{S{u))), 

which yields the claim. □ 
Lemma 4.4. For every e > 0, there exists M(e, C) > such that 

ll^d4m')<e\\VMl + M\\{i;,u)\\l (4.27) 

for any u G UngnXn, lo = X/ x u, and ip = X/ x uj = —Au. 
Proof. Recall (cf. (glD) that 

i^^lVp) = -vr.((V;,i^)V'ii)-^(V'",V'")-^IV'^l' 

+2(V",Vr(^,zy)) + (V" X {V xij)\\,i^). 
Integrating the equation over and using the boundary data in Lemma 14.31 yield 

+iy"(V",5(u))-2^7r(V'll,5(n)) 
2 

~c 

The first integral vanishes by Stokes' theorem applying to the surface F. Using the Holder 
inequality, it follows that 

llMli^f) < M\\i,m,^^^{l + \\S{u)U^r) + \\{V\Vu,u)h2(^r))- (4.28) 

All the boundary integrals on the right-hand side can be estimated via the Sobolev imbed- 
ding: 

\\{V\,Vu,u)\\L2^r)<e\\VM2 + M\\{i;,u)\\2 

and 

||5(n)||i2(r) < e||V5(n)||2+M||5(n)||2<e||VxV'||2 + M||(^,^x)||2, 



^(t/;", (V^ • n)) + 4^(V^II, (V^ • niu))). 
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where the second inequahty follows from (j4.16p . Again, by the Sobolev imbedding, 

Plugging these estimates into (|4.28|) . we conclude 

llduW) < (^llVVlb + C^ll^lb) {e\\V\h + C\\i^,u)\\2 + l), 

and (OTj) follows. □ 
Theorem 4.4. Let u G UNef>i^N, uJ = V x u, and tp = — Au. Then 

||V\||2 < M(||V^||2 + \\u\\h2) 

so that 

||n||H3 <M||(VV^,V,tx)||2, (4.29) 

where M > is a constant depending only on ( and the domain 0, which may be different 
in each occurrence. 

Proof. First, we need to apply integration by parts to reduce the L^-norm of the third 
total derivative V^u into the L^-norm of the total derivative Vtp, plus some boundary 
integrals, which can be dominated by the use of boundary conditions on u and the facts 
that V • u = and V • V' = 0. Let us carry out this step in the ordinary coordinate system, 
and thus dj = are the coordinate differentiations. Then 

iv^ni^ = {d^djdku^){didjdku^), 

where, as usual, the repeated indices are summed up from 1 to 3. Therefore, by using 
integration by parts twice, 

l|v'u||i = / {dk^'){dk^')+ [ {dk^'){djdku'){dj,ij)+ [ {djdku'){didjdku'){di,iy) 

JQ Jr Jr 

IVV'I' + ^ {dki^') {d.dku') {d,,u) + i ^ a,(|V\|2) , (4.30) 

where ■0' = — Au' has been used. We first handle the last boundary integral | Jp (| V^np). 
To this end, we use a moving frame (Vi, V2, V3) so that V3 coincides with the unit normal 
u and Vj, i = 1,2, are (local) tangent vector fields, so that we may perform integration 
by parts on the surface F for Vj, i = 1, 2. Then 

J=IM\^^^\^) ~ (V,Vfc^,')(v,VfcV3^x'), 

where ~ means that the difference between the two sides contains only the quadratic terms 
involving u and its 1st and 2nd order derivatives, the second fundamental form vr and its 
derivatives. These terms come from the exchange of orders by applying the derivatives 
Vj. According to the Ricci identity, the difference after the exchange of VjVj to VjVi is 
a term which has order 0. More precisely, 

3 

v,v,r'- v,v,r' = J^c^.T', 

1=1 
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where T is a vector field and 1 can be a multi-index, depending on the order of the vector 
field T. The boundary integrals of these lower-order terms can be controlled by 

C7||(VVVn,n)||i.(r), 

which, in turn, can be dominated by 

e\\VMl + M\\{'^,n)\\l. 

Note that, in this step, no boundary condition on u is required. With this principle, 
without using any boundary condition on u, we may re-group [VjVkU^^ (VjVfcVau') 
into the following items: (VjVfcU^) (Vj-V^Vsu^) , (V^Vsu^) (VfcViw^), (Vlw^) (V|n^), 
(VjVfcu') (VjVfcVsu'), (VjVsii') {VjVlu^), and (V|ii') (V^u'), where ah the other re- 
peated indices run through from 1 to 2. Indeed, 

j,k=l,2 l,j,k=l,2 
l,j=l,2 k=l,2 

-f(Vi«3)(V3Vin3)+5;(Vin')(Vin'). 

1=1,2 

To estimate the boundary integrals of the right-hand side, we use the kinematic condition 
(jl.3p . the divergence-free condition: 

Vsn^ = -Yl ^^^^ ^sV'' = -Yl (^-^l) 

a=l,2 a=l,2 

and the definition of tp: 

Vin" = Au'-Y: V^u" = -V" - E (4-32) 

a=l,2 a=l,2 

Both ()4.3ip and ()4.32h hold on Q. By using these relations, we can rewrite J as follows: 

+ {2ViVau'' + V^n' + (V,V') + (Vgn' + (VsV' - V^^^) 
+V'v2(V3n' - Viu^) + E (V.Vfcu') {V,Vk{Vsu' - Vm^)) 

l,j,k=l,2 

+ E (V^0V^(V3n'-V,n3), 

«,fe,a=l,2 

where the other repeated indices are added up through 1 to 2. The terms in the square 
brackets are to be replaced via the corresponding boundary conditions, so that the orders 
of taking derivatives for these terms are reduced by 1. Therefore, all the terms, except the 
first two, are quadratic forms of u and its 1st and 2nd order derivatives, and those of vr, 
so that these terms are dominated by 

C\{V^u,Vu,u)\^. 

For the first two terms, we perform integration by parts on T. The second boundary 
integral on the right-hand side of (j4.30p can be estimated similarly. 
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Let I = Jp9iy(|V^«p). Integrating the above equation and using integration by parts 
on the surface F yield 

/ < ^ {2^'^ + 3Vfc V„tx« + Vlu'') Vfe V' + / (V'^ + V&^O (VaV-' - V/ V') 

+e\\V\\\l + M\\{i;,u)\\l 

Therefore, using the boundary conditions for = l(ip,v), V ^il)^ — V iil)"^ (which is (VX'^)II), 
Va^' — Viu^ (which is (V x n)"), together with the Sobolev imbedding, we have 

I <e\\V^u\\l + M{eX)\\{^,u)\\l 

□ 

Corollary 4.2. There exists M{Q > depending only on ( and such that 

M\\{VtP,tP,u)\\2 < \\u\\h3 < M-i||(VV',V',^)l|2 (4.33) 
for any u G UNeN^N, where ip = — An. 

4.4. The Stokes semigroup. The self-adjoint operator S on K2{^) has a spectral de- 
composition 

S= r XdEx, 

J — oo 

where {Ex : A < A} is the left-continuous family of the projection operator Ex on the 
space spanned by {ofe : > A}. According to the spectral theory of self-adjoint operators, 
u G K2{^) belongs to -Df(S') if and only if 



r X^d{Exu,u) = f2xl( [ {ak,u)) 



2 

< OO, 



and D{£ + M)=D {^-S + A/) = K2{n)nH^{n); and u G ^2(0) is in D{S) if and only 
if 

r-A 



/ Xd{Exu,u) = y2^k( [ {ak,u)) < 00. 

J -00 ZTn ' 



k=0 

In this case, 

S{u, = - ^ Afc ( / (afe, «)) . (4.34) 
fe=o 

We are going to show the following estimate that plays an important role in the proof 
of the existence of strong solutions to the Navier-Stokes equations. 

Theorem 4.5. For any £ > 0, there exists M{£, Q such that 

\\VxPn{u)\\1<{1 + £)S{Poo{u),Poo{u)) + M\\u\\1 (4.35) 

for any u G Lp'{Vt) and any integer N . 
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Proof. Recall that 

N „ N 



PNiu) = ^ak / {ak,u) = ^ak / (afc,PooN> 

7 — n 7 — n ^ 



A:=0 " A:=0 

SO that PAr(n) G Z:'o,c('S'). According to ([Q]) . 

||VxP^(tx)||2 = ||VPiv(n)i + y"7r(P^(tx),P^(tx)) 

< (l+e)||VP^(n)||i + M(.,C)||P7v(^)i 

< (l + e)£:(Pjv(n),P7v(t^)) + M||n||^ 

= -(1 + e) [ {S {Pn{u)) , P,v(n)) + M\\u\\l 



where the second inequality follows from (|4.13p . However, 



AT 

S'(PAr(n)) = ^ AfcOfc / {ak,u). 

k=0 ''^ 



Denote integer A^'o > such that > \q > ■ ■ ■ > Xno >0> \nq+i > • • • • Then we find 
that, when > Aq, 



[ {S{Pn{u)),Pn{u)) = -y2>'k [ {ak,PN{n)) [ {ok, 
Jn Jn Jn 

= -'^^ki {ak,Poo{u))) 

1 n ^ 



k=0 



< 



y]Afc( / {ak,Pooiu)))' 



k=0 

£{Poo{u),Poo{u)), 



and, while A^ < A'o - 1, 



/ {S{PNiu)),PN{u)) = -Y^Xki [ {ak,Pooiu))f 

oo „ TVo „ 

^ -y.M {ak,Poo{u))Y + V Afe( / {ak,Poo{u))y 



k=0 " " fc=Af+l 

= £{PUu),PM)+M\\u\\l. 

This arrives the result expected. □ 

Corollary 4.3. For any e > 0, there exists M{e, C) > depending only on e, Q, and Q 
such that 

||(V X Pn{u),VPn{u))\\1 < M||(V X u,u)\\l (4.36) 
for any u G H'^{Q) and any integer N . 
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Proof. Note that Pooiu) € K2i^) n H^{n) so that 

V X Poo{u) = V X Poo{u) = V X n 
and the estimate follows from 

£:(Poc(n),Poo(n)) < M\\{VPo,{u),u)\\l 
which yields (j4.36p . Similarly, the estimate for VP/v(ii) follows from (j4.36p and (|3.4p . □ 

5. Existence of Weak and Strong Solutions 

In this section, we consider the initial-boundary value problem (jl.ip - (jl.4p . where C > 
is a constant. The minimal requirement on the initial data is that uq £ (^)- 

First we introduce the notion of weak solutions to the initial-boundary problem (jl.ip - 
(jl.4p . We say that a vector field u{t,x) is a weak solution of (jl.ip - (jl.4p with slip length 
C, provided that u{t, x) satisfies the following conditions: 

(i) For each t > 0, u{t, •) G K2{n), and u G L'^{[0,T]; H^{n)) for any T > 0; 

(ii) For any smooth vector field (p(t, x) with ip{t, •) G K2{^1), 

(n(r,.),</^(r,-)> 

n 

{uo,fo)+[ [ {u{t,-),dtf{t,-)) 
Jo Jn 

(V X n, (n X + /iV X (/?))- ^ / {u,ip) + 2fi TT{u,ip); (5.1) 



Jn 



CJo 



(iii) The energy inequality: 



||n(T,-)||i + 2/u^ \\Vug + 2fiJ^ ^(i|n|2-7r(n,n)) < llnolli. (5.2) 

Equation (j5.ip is obtained by integrating (jl.ip and performing formally integration by 
parts. 

5.1. Construction of global weak solutions. Notice that (5, D^(5)) has the eigenval- 
ues A > Ao > Ai > • • • > A,i — > — oo and the eigenvector functions {a„ : n = 1, 2, • • • }, 
which form an orthonormal basis of i^2(^)- Then San = ^nO-n, and solves the Stokes 
equation: 

Aa„ - Vpn = Xnttn, V • On = 0, (5.3) 

subject to the kinematic condition (jl.3p and the Navier's C-condition ()1.4p . 
Let 

oo „ 

u{t,-) =y2ck{t)ak with Ck{t) = / (afc,u), 
k=i 

be a solution of the Navier-Stokes equation (II. ip with initial data uq- Multiplying by 
and integrating over Q, we obtain 

dtCk = M / {ak,^u) - / (afc,^i • Vn) = /uAfcCfc - V] QCj / (aA,.,aj • Voj). 
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Thus, for each integer N, we solve the Cauchy problem for the system of differential 
equations: 

N 



dt 



Cfe = /xAjfcCfe - CiCj / {au.ai ■ Vaj) (5.4) 
i,j=o 

Cfe|t=o = / {ak,uo)- (5.5) 
Jn 



Define 

N 

.N 



{t, •) = ^Cfc(i)afc. 



u 

k=0 

Then u^{t, •) G Do,c('S') for f > and satisfies the evolution equation: 

AT 



dtu'' = ^iS{u'') -Y^aJ {ak,u'' ■ Vu^). (5.6) 

k=i ■'^ 



Therefore, V • = 0, {u^)\ = 0, and (V x u^)^^ = -i(*u^) + 2( * 7r(n^)) for t > 
Now we make the energy estimates for and . First, it is easy to see that 

TV , AT 



fe=l k=l 



Since 

AT 



k=l 

we have 

dt 



^||n-|l2 = -2^/ \Vu^f- /t/^-V(|n^|2)-^ / |n^|2 + 2/. /7r(./', u^^). 

Jn C Jr Jr 



Since V ■u'^ = and (■u'^)-'- = 0, which implies ■ V{\u^\'^) = 0, we obtain the energy 

balance identity: 

^llu^lli + 2;x||V..^||i = -^||u^||i.(r) + 2m j^i^{u\u^). (5.7) 



Therefore, we have 



u^{T,.)\\l + 2^^|^ \\Vu^\\l + 2nj^ j^{h,u^\^ - ■,{u\u^)) = \\u^ \\l < \\u4l (5.8) 
Using the Sobolev embedding inequality: 

j^\u\^ <e\\Vu\\l + C{e)\\u\\l 

|^7^(^^u^)<i||V^^||i + C(, 



we have 



e)\\u^\\l- 
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Then, from (jS.Sp . we have 

lk^(r,-)lli + ^ri|Vn^i + ^ r [ \u^\'<\\u4l + C r\\u^is,.)\\l (5.9) 
Jo s Jo Jt Jo 

The Gronwall inequahty and (|5.9|) imply that •)||2 and /g ||Vm^||| are uniformly 

bounded in t, C, and A^. The apriori estimate (jS.Sp also ensures that, for each integer A^, 
system ()5.4p has a unique solution for all t > 0. Then we conclude 

Theorem 5.1. Let uq G K2{0,). Then, for any T > 0, the family {u^{t,x)},0 <t <T, 
is weakly compact in the space L'^{[0,T]; K2{0,)) so that it has a convergent subsequence 
that converges to a vector field u € L'^{[0,T]; K2{^)), and the limit function u{t,x) is a 
weak solution to problem (ll.l[ )- (ll.4l) . 



Furthermore, we have 

\WdT.-)\\l + l^f I l|Vnc(t,-)lli + 7^ f I Kl'<IKIli, (5.10) 
Jo Jn ^ Jo Jr 

where we have used to indicate the dependence on the slip length C > 0. The uniform 
energy estimate (|5.10p implies that the family U(^{t,x) is pre-compact in L'^{[0,T]; K2{0,)) 
and hence there exists a convergent subsequence (still denoted) u, when C ^ 0. 

From (j5.10p . we have 



|2 

l25 



which leads to 

rT 



Jo Ir 



= 0. 



'r 

This implies that the limit u{t,-) is subject to the no-slip condition for almost all time t. 
On the other hand, when — > oo, we obtain that there also exists a subsequence (still 
denoted) U(^{t,x) converging to u{t,x) in L^([0, T]; i^2(^^)) such that u{t,x) is a solution 
to (II. ip subject to the complete slip boundary condition: 



in the weak sense. 



2(*7r(n)) (5.11) 



Theorem 5.2. Let U(^{t,x), < t < T, be a weak solution to problem (ll.ip - p.4p con- 
structed in Theorem \5.1[ Then 

(i) when C ^ 0, there exists a subsequence (still denoted) U(^{t,x) converging to u{t,x) 
in L'^ {[0,T]] K2{0.)) such that u{t,x) is a solution to (II. ip subject to the no-slip condition 
for almost all time t; 

(ii) when C, — > oo, there exists a subsequence (still denoted) Ui^{t,x) converging to u{t,x) 
in L^([0, T]; i^2(^)) such that u{t,x) is a solution to (11. ip subject to the complete slip 
boundary condition (15. lip . 

The nonhomogeneous vorticity boundary problem related to (15. lip has been investigated 
in [7]. 
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5.2. Strong solutions. In this section, we prove that there exists a strong solution to 
problem (jl.ip - (jl.4p for small time. To this end, we develop the L^-estimates for up to 
second-order derivatives, uniformly in N. More precisely, we prove the following: 

Theorem 5.3. Let uq G K2{^) n H'^{9). Then there exist T* > and M > depending 
only on Q, e, fi, Q, and \\uQ\\ff2 (but independent of N) such that 

\\u^it,-)fH2 + \\dtu^it,-)f2<M. (5.12) 

Proof. For simplicity, we write u = given by (|5.4p . Let a; = Vxn and i/; = Vxtj = — Au 
as usual. Recall that u fulfils the evolution equation: 

N 



dtu = fiS{u) -S^CLk / {ak,u-Vt 



where, for t > 0, u{t, •) G Dq^(^{S). Taking the t-derivative yields the evolution equation: 
dtut = nS{ut) -^Ok / {ak,ut -Vu) - ^Ok / {ak,u-Vut), (5.13) 

k=l k=l 

and ut{t,-) G Dq^(^{S). Therefore, V ■ ut = and ut again satisfies the same boundary 
conditions as those of u. Using the evolution equation (j5.13p . we have 

^ll-Uills = 2fi f {S{ut),ut) - 2 f {ut,ut -Vu) - 2 f {ut,u-Vut) 



dr 

= 2ii I {S{ut),ut) -2 I {uuUfVu). (5.14) 



Integration by parts yields 

(Ant,ut) = -\\Vut\\l+ I {utx{V xut),v) + I {ufVuuv) 



-ll^^tlli j {utx {*ut),v) + 2 J {utx {*TT{ut)),I^) - j 1T{ut,Ut) 



Substitution this into (|5.14p leads to 

^\\ut\\l = -2fi\\Vut\\l ~ ■ ~ y / 1^*1^ + '^'^ J^^i^t, Ut). (5.15) 



dt 

It follows that 

d II,. i|2 < 2^ 



^Il«tll2 S --WutU.^r) 



2/i||Vut||2 + 411^(11211X77^1100 + 2;U j '^iut,ut) 



< IKtlli2(r) - /"llVntlli + e\\V\f + C {\\ii;,u,ut)f2 + \\ut\\l) , (5.16) 

where the Sobolev's imbedding has been used and C > is a constant depending only on 
the domain Q. 

Next we deal with An. The evolution equation for u{t, •) may be written as 

dtu = fiS{u) — Pn {u ■ Vu) . 
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Together with the vector identity u ■ Vu = ^V|np — u x uo, we have 

dtu = fiS{u) + PNiuxuj). (5.17) 
Since V x S{u) = V x (An), taking curl (twice) to both sides of equation (j5.17p yields 

dtuj = IJL/S.UJ + 'V X Pn{u X Lo), (5.18) 

and 

dt'ijj = nAip + V xV X PNiuxuj). (5.19) 
It follows from (f5T9]) that 

^\ml = 2fi [ {AiP,iP)+2 [ {V xV X PNiuXLo),^P). (5.20) 



Integration by parts leads to 
and 

(V X V X Pat (n X w) , -0) 



-2^||VV'||i + /i^9, 



{V X Pn (u X Lo) ,V X ip) - J X {V X Pn{ux uj)) , u) 

{V X Pn {u X uj) ,V X ip) - I iiJ X { - ]:{*Pn (u x lo)) + 2(*7r {Pn (u x uj)))),i^) 

Jr C 

(V X Pat (u X w) , V X V') + 7 / {ip, Pn {u x t^)) - 2 / vr {ip, Pn {u x u)) . 
n C Jr Jr 

Therefore, we obtain 



^llVlli = -2^11 VVlli + 2 ^(V X Pat (nxc^),Vx^)+/.^ 5, (|^|2) 

+^ j^{i;,PNiuxu;)) -4j^Tr{i;,PNiuXio)). (5.21) 
Using the Holder inequality, one obtains 

^llV'lli < -2^||VVi + 2||V X P^(n XL^)||2||V X ^||2 



+fi / d,{\ip\^) + MU\\L2^r)\\PN {u X oj)\\L2^r)- (5.22) 



d,{\^f) <e\\V\g + M\\{^,u)\\l 



The first boundary integral /p<9,y(|^p) can be estimated by using Lemma 14.41 to obtain 

1 
2 

The product of last two boundary integrals in ()5.22p can be estimated via the Sobolev 
imbedding to yield 

^IIV'llL2(r) \\Pn {u X to) ||i2(r) < {eWViPh + MW^Ph) (e||VPjv (n x a;) h + M\\u x coh). 
Plugging these estimates into (j5.22p , using Corollary 14.31 and the estimate 

\\uxuj\\2 < M||u||^i = M||(V',n)||2, 



I < -2^l\\y^l + e\ml,+e\\V^u\\l + M\\VPN{u^u:)h\mH^ 



26 GUI-QIANG CHEN AND ZHONGMIN QIAN 

and rearranging the inequality, we obtain 
d_ 
dt 

+M(||(V^,^)||i + ||(^,n)||^). (5.23) 
Finally, we use Corollary 14.31 to obtain 

\\V Pn {u X LO) \\2 < M \\V X {u X Uj) \\2 + M\\u X Ll!\\2 

< M\\{uj ■ Vu, u ■ Vlo, u X u;)\\2 < M\\u\\jj2- 

Then we conclude 

j^ml < -^^llV^lli + M{\\{i;,u)f2 + ||(^,n)||^) (5.24) 

for some constant M depending only on C, /i, and 0. 

Let F = \\{ip,u,ut)\\2- Combining (j5.7p and (j5.16p with (j5.24p . we obtain the following 
differential inequality: 

j^F < -^11 V(^, n, ut)\\l - y II {u, ut)\\l2^r) + MiF + MaF^ (5.25) 
for some constants Mi and M2 depending only on C and 0,. In particular, we have 

^F < MiF + M2F^. (5.26) 

Since 

Wuth < /"l|5'(n)||2 + ||PAf {u X 1^) II2 < 2/i||An||2 + ||n x t<j||2 < 2/i||An||2 + ||n||2||c(;||2, 
then 

F(0) <C(/u,0)||txo||^2 
for some constant C (/i, f^) depending only on /j, and 17. 

Let p be the solution on [0,T*) to the ordinary differential equation: 

p' = Mip + M2p\ piO) = C{p,n)\\uofH2, (5.27) 

where T* > is the blowup time of p. 

Then inequality (j5.26p together with the fact that -F(O) < p{0) implies that F{t) < pit) 
on [0, T*). This completes the proof of Theorem 15.31 □ 

Theorem 5.4. Let uq £ K2{^) ^ H'^{Q). Then there exists T* > depending only on 
C,p, ^, o-nd \\uq\\jj2 such that there is a strong solution u{t,x) of the initial-boundary value 
problem ([TT]) - (frt]) up toT* >0. 

6. Inviscid limit as p ^ 

In this section, we analyze the inviscid limit of the solutions u^{t,x) of the initial- 
boundary value problem (|l.ip - (|1.4p . 

Let u{t, x) be the unique smooth solution of the initial-boundary value problem of the 
Euler equations: 

dtu + u ■ Vu = —Vp , 

uiO,-)=u{0), ^^-^^ 
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up to time T* > 0. Notice that all solutions to problem (jl.ip ~ (jl.4p . fi £ (0, ^uq] for 
some fiQ > 0, subject to the same boundary conditions: = and 

(V X = + 2(*^(n'^)), 

while the solution u of (j6.ip satisfies only the kinematic boundary condition and is inde- 
pendent of the viscosity constant fi. 

Theorem 6.1. Suppose that, for all fi £ (0, ^uq], a unique strong solution of the initial- 
boundary value problem (|l.ip - (|1.4p and a unique strong solution u G if^(Q) to the initial- 
boundary value problem (j6.ip ) both exist up to time T* > 0. Then there exists C = 
C{iJ,o,T, ||ii||L2([o,r];_ff2nVKi'°=(Q))? independent of fi, such that, for any T G [0, T*], 

sup •) Olb < C/i ^ as /iiO, (6.2) 

0<t<T 

and 

rT 

\\V{uf-u) {s,-)gds < C. 

It follows that the whole solution sequence u'^ of (|l.ip - (jl.4p converges to the unique solu- 
tion u{t, x) of the initial-boundary value problem (|6.ip in as u ^ 0. 

Proof. Let = — u. Then satisfies the following equations: 

^tv^' = fiAvf" - (vf" + u) ■ Vv^' - VP'' -v^" ■Vu + fiAu, , . 

V■v^' = 0, ^^--^^ 

and the initial condition: 

vf'{0,-)=0, (6.4) 

where = p^ — p. Since both and u satisfy the kinematic condition (|1.3p . so does v^. 
Thus, by means of the energy method, we obtain 

d 



2fi [ {v^'^Av^')- [ {v^ + u,V{\v^'\^)) -2 [ {VP^',v>') 
Jn Jn Jn 

-2 [ {vf" ■Vu,vf') + 2n [ {Au,vf'). 
Jn Jn 



Integration by parts in the first three integrals leads to 



d 



{v'' ■Vu,v'') + 2n / (An,?;'') 
Jn 

2 [ {v'' ■Vu,v'') + 2n [ {Au,vf'), 
Jn Jn 



where b = -^{* {v'' + u)) + 2( * 7r{v^' + u)) - (V x n)". 
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Furthermore, we use the fohowing estimate: 

{v^ xb,iy) < \\b\\L2(^Y')hf,\\L^r) < C(lhllHi(r) + Il^''lli2(r)) 



T 

to obtain 

+2||Vn||oo|b''||2 + Ml^uhlKh- (6.5) 
Finally, we use the Sobolev imbeddings: 

ClKWli^j^) < \\VVf,\\l + C\\vjl, ll'ull^i(r) < C\\u\\jj2(^^) 
to establish the differential inequality: 

^ll^/^lli + < C (||Vn||oo + Ho) \\v^\\l + /^||'u||^2(s^). (6.6) 

Gronwall's inequality implies that 

\K{t,-)\\l < /u/o*e^(/.'ll^«(-'-)ll-'^-+^o(*-^))||n(s,.)||^2(o)ds =: Cfx, (6.7) 

where C > depends only on fj,o,T*, and \\u\\L'^[[o^T];H^nW'^'°°{Q)- Using this and (16. 6p . we 
further have 

ft 

\VVf,{s,-)\\lds < Ch. 

'0 

This completes the proof. □ 

In order to ensure the convergence of to u in the strong sense (say, i?^(n)), a 
necessary condition is that u must match the Navier's ^-condition (jl.4p . 
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